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PREFACE 

This  report  elaborates  and  extends  some  of  resiilts  of  the  Theory 
of  Cusped  Geometries,  I.  General  Survey,  by  Harold  Grad  [i].]  .      The 
computation  of  particle  losses  through  a  line  cusp  is  performed  under 
the  assumption  that  there  is  a  good  approximation  to  a  sharp  boundary 
separating  the  field-free  plasma  from  the  vacuum  magnetic  field.   A 
rough  calculation  along  the  lines  of  this  report  was  first  presented 
at  Sherwood  meetings  at  Los  Alamos  in  June  1955  (no  report  issued), 
and  at  Princeton  in  October  1955,  [2] , 

Drc  E.  Rubin  was  very  helpful  in  the  preparation  of  early  versions 
of  this  reporto 
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1.   Introduction 

The  pioppose  of  this  report  is  to  compute  the  losses  of 
particles  through  the  cusps  from  a  cusped-shaped  fully  ionized 
plasma  separated  from  a  vacuum  magnetic  field  by  a  thin  boundary 
layer.   (See  Figure  1;  the  thin  boundary  layer  is  not  shown.) 
To  find  the  losses  one  must  use  a  more  realistic  model  than  that 
of  the  simplified  magnetohydrodynamic  picture  of  a  sharp  boundary 
between  plasma  and  vacuum.   Among  other  matters,  the  equilibrium 
and  stability  of  cusped  geometries  is  treated  in  [5,6,?]. 


Figure  1 
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We  first  consider  a  two-dimensional  problem  concerning  the 
motion  of  particles  in  cusped  domains  satisfying  a  certain  unconven- 
tional reflection  law,  "soft-cushion  billiard  ball"  reflection,  which 
simulates  the  finite  thickness  boundary  layer , cf .  [2] .  We  determine  the  loss 
rate  of  particles  from  the  cusp  according  to  this  model,  and  finally 
we  demonstrate  the  connection  between  this  model  and  the  actual  motion 
of  charged  particles  in  the  plasma. 
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2„   Orbits  in  the  plasma 

To  fix  our  notation,  see  Pig.  2,  we  take  a  rectangular  x,y  coordinate 
system  in  vjhich  the  axis  of  the  cusp  is  the  x-axis,  with  the  cusp 
at  the  origin;  the  domain  of  Interest  extends  sjnnme  trie  ally  about 
the  negative  x-axis „   We  denote  by  \|/  the  angle  between  the  positive 
X-axis  and  the  direction  of  motion  of  a  particle,,   When  a  particle 
crosses  the  cusp  axis,  if  it  is  advancing  toward  the  cusp,  we  have 
I  ij;  I  <  -^  J  whereas  if  it  is  moving  away,  |  \|;  |  >  S- « 

We  make  the  following  general  assumptions: 

In  the  interior  of  the  domain  a  particle  moves  on  a  straight 
line,   VJhen  it  flies  out  of  the  domain  into  the  boundary  layer,  it  is 
returned  shortly  via  some  trajectory  which  need  not  be  specified  in 
detail  now.   Then  it  continues  inside  the  domain  on  a  straight  line 
in  the  direction  it  had  at  the  instant  of  its  return  until  it  emerges 
again,  is  returned  again,  and  so  forth.   Our  analysis  is  made  close 
to  the  tip  of  the  cusp,  where  the  boundary  curves  of  the  cusp  are 
very  nearly  parallel  straight  lines. 

Upon  reentry  into  the  plasma  domain  a  particle  moves  in  that 
direction  which  it  would  have  if  it  had  undergone  specular  reflection 
at  the  magnetic  wall, 

A  particle  in  two  successive  flights  outside  of  the  cusp  domain 
advances  (in  the  x-dlrection)  a  distance  r,  +   r^   =   2r(\lf);  r(\|/)  is 
called  the  range,  (see  Figure  2b),   The  range  depends  continuously 
on  "if    in  such  a  way  that: 

(a)  r(\l/)  >  0  0  <  \1/  <  I  , 

(b)  r(7i:-\l/)  =  -r(\|;)      0  <  ^  <  it  , 

(c)  r(-\l/)  =  r(\lf)        0  <  \lf  <  Ti  . 
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Pigiire   2a 


to  cusp 


Figure   2b 
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This  prescription  completely  determines  the  motion  of  a  particle, 
provided  the  particle  path  never  extends  beyond  x  =  0,   When  the  latter 
situation  arises,  we  say  the  particle  is  lost  through  the  cuspo   Thus 
■we  regard  as  lost  all  particles  not  tiorned  back  before  reaching  the 
tip  of  the  cusp.   What  happens  to  particles  which  are  lost--whether 
they  return  to  the  interior  of  the  plasma,  or  to  the  boundary  layer-- 
Is  a  complicated  problem,  which  Is  not  considered  here.  However,  see  [I4.]  , 

In  the  vicinity  of  the  cusp,  the  boundary  curves  of  the  cusp  are 
nearly  parallel  straight  lines,  and  the  range  r(\l/)  is  small  compared 
to  the  dimensions  of  the  cusp  domain.   If  we  denote  by  /\\|/  the  change 
in  \j/  after  two  successive  reflections  from  the  exterior,  then  it 
follows  that  lA^^I  «  \^\'-      Thus  we  have  to  good  approximation 
(cf.  Figure  2b),  for  ^   positive, 

^x  =  l\.j   cot  \|/  +  2r(\lf) 
and 


and  we  find 


Ax 


-^     =   -y  cot  \|/  -  i  r(\l/)  . 

At  ^ 

For  small  A^s  "®  ^^^   replace  this  difference  equation  by  the  corres- 
ponding differential  equation" 

^  =   -y  cot  \|;  -  I  r(\lf), 

which  can  be  integrated  immediately  to  yield 

"This  analysis  is  approximately  equivalent  to  using  an  adiabatic 
invariant  which  generalizes  the  magnetic  moment;  cf,  [3]  and  the 
appendix  to  [5] » 
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Tl/2 


y  = 


sin   \1/ 


I  r(\[')    sin   \|/   dij/   +   y(|) 


We  sketch  some  of  the  integral  curves  in  Figure  3< 


•-y 


Figure  3 


Note  that  there  is  a  critical  integral  curve  for  which 
y(-5-)  =  0.   We  denote  this  curve  by  \|/  =  \l'"(y);  the  inverse  of 


this  function  is  given  by 


Tl/2 


y  = 


2  s 


^^   ^-;:-   I    r(\l/)  sin  \|;  d\i;  . 
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This  critical  curve  \|;  =  \|/"(y)  has  the  property  that,  for  a  given 

value  of  J,   (the  half -width  of  the  cusp),  all  particles  starting  with 

i^l  -^  H'"(y)  will  be  returned  before  reaching  the  cusp,  whereas  those 

starting  with  |  \|/ ]  <  \}f"(y)  .^fill  be  lost  through  the  cusp^   We  remark, 

for  future  use,  that  2y\lf'"  tends  to  a  constant  as  \i/"    tends  to  zero, 

and  y  tends  to  oo ;  we  denote  this  constant  value  by  r'"',  where 

n/2 
(1)      r'"  =     r(\(;)  sin  \t(  d\|;. 


o 


Incidentally,  in  the  special  case  of  specular  reflection,  i.e 
r(\(f)  =  0,  the  integral  curves  are 


7   =  y(f  )/sin  \1/  , 


and  the  critical  integral  curve  is  y  =  0,   Thus  in  this  case  our 
theory  yields  the  (correct)  result  that  the  only  particles  lost  are 
those  aimed  exactly  along  the  cusp  axis  towards  the  tip* 

It  is  important  to  note  that  our  conclusions  are  independent  of 
the  shape  of  the  cusp „  The  equation  y  =  y(x)  of  the  plasma  boundary 
does  not  matter;  one  only  needs  to  know  the  range  r(\|;). 

Having  determined  the  path  of  a  particle  moving  in  a  plane  cusp- 
shaped  domain  in  accordance  with  the  laws  of  motion  specified  above, 
we  now  consider  a  general  three-dimensional  motion  in  a  cusp-shaped 
wedge  in  x,yj,  z-space.   Now  Figiire  2  represents  a  section  of  the  wedge 
for  z  constant,   (Note  that  we  are  now  dealing  with  a  line  cusp  J  We  sup- 
pose  that  a  gas  consisting  of  particles  is  contained  in  the  cusp-shaped 
domain.   We  shall  compute  the  rate  of  loss  of  particles  through  the 
cusp. 
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In  the  physical  model  described  below  it  will  be  seen  that  the 
motion  of  a  particle  is  such  that  the  projection  of  its  trajectory 
on  the  x,y-plane  satisfies  the  general  assumptions  made  above.   The 
components  of  velocity  in  the  x,y,z  directions  are  denoted  by  u,v,w 
respectively.   We  also  write  u  =  U  cos  ^   and  v  =  U  sin  \|;.   In  the 
Interior  of  the  cusp-shaped  domain  we  assume  that  the  speeds  U  and  w 
remain  constant  for  a  particular  particle  in  the  course  of  its  motion,, 
The  range  r  now  depends  on  U  and  w  as  well  as  on  \|/;  ioe.  r  =  r(\|f,U,w), 
We  shall  ass\ime  that  the  distribution  of  particle  velocities  is  isotropic 
far  enough  from  the  cusp. 

We  remark  first  that  the  distribution  of  particle  velocities 
does  not  remain  isotropic  as  the  cusp  is  approached.   If  we  calculate 
the  rate  of  flow  through  a  section  of  width  2y  at  some  given  x,  of 
those  particles  which  have  |  \|f  |  <  ^'^ij),    we  find  that  this  quantity  is 
independent  of  y.   Hence  it  is  easiest  to  perform  the  calculation  at 
a  relatively  large  value  of  y  where  the  velocity  distribution  is 
isotropic.   This  loss  rate  will  now  be  computed  in  terms  of  r(\);,U,w) 
or  r'"'(U,w)  o 

The  isotropic  velocity  distribution  function  f(u,v,w)  is  given. 
The  flow  rate  per  unit  area  in  the  positive  x-direction  is 


u  f   du  dv  dw     =      j      u  f  U 

u>o  kl<f 


dU  d\|/   dw 


9 


and  the   loss  rate   per  unit   area,    determined  by    |  \li  |    <  \|f'"',    is   easily 
seen   to  be,    for  ^'"    small, 


2\l("(U,w)   U'^f  dU   dw. 
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Therefore  the  loss  rate  per  unit   length  of  line  cusp  is  given  by 

Q  =  2    r'^U  f  U  dU  dw 
or 
(2)    Q  =   1   |rvl  f  du  dv  dw  , 

u  >  0  -  , 

Within  a  numerical  factor  of  order  xinity,  the  mean  value  of  r' 
represents  the  effective  size  of  the  "hole"  at  the  cuspo 
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3,   Range  f-unction 

We  shall  now  calculate  the  range  r,  which  was  found  to  be  basic 
in  computing  particle  losses  through  cusps.   To  do  this  we  need  to 
choose  some  model  for  the  boundary  layer  between  the  plasma  and  the 
vacuum.   We  adopt  the  simple  model  of  a  magnetic  field  B,  which  jumps 
discontinuously  from  zero  to  its  vacuum  value  together  with  a  constant 
electric  field  E,  caused  by  charge  separation.   We  can  either  compute 
a  roughly  self-consistent  value  of  E,  e.g.  by  equating  the  average 
penetration  of  ions  and  electrons  into  the  bo\indary  layer,  or  else 
leave  open  the  value  of  E.   The  latter  coiorse  is  appropriate  in  case 
there  are  stray  electron  orbits  near  the  plasma,  but  not  entering  it»" 
The  number  of  such  external  electrons,  which  influence  the  value  of  E, 
would  have  to  be  determined  by  an  analysis  of  the  whole  system.   Order 
of  magnitude  estimates  have  been  made,  cf,  [l\.]  . 

In  the  boundary  layer  we  introduce  (x,y,z)  coordinates  as  Indicated 
In  Figure  l\..      The  x-axis  is  directed  towards  the  cusp,  and  the  y-axls 
is  in  the  outward  normal  direction  to  the  plasma.   Thus  we  use  a       i 
different  coordinate  system  in  the  other  boiindary  layer.   The  components 
of  velocity  are  (u,v,w),  as  before.   It  is  convenient  to  Introduce 


V 


=  /u+v+w  and  v'    =   e/B    . 


Now  suppose  a  particle  enters  the  boundary  layer  atx=y=x=0 
with  initial  velocity  (u,v,w).   It  is  a  simple  matter  to  find  the 


'cf.  M.  Rosenbluth  [1],  who  gives  a  more  general  self-consistent 
analysis,  but  treats  a  special  velocity  distribution.   However,  a 
more  precise  self-consistent  treatment  of  the  boundary  layer  is  not 
warranted  here  because  of  the  uncertainty  as  to  the  exact  value  of  E. 
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E  =  0 
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boundary  layer 
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trajectory  of  a  particle  with  mass  m  and.  charge  e  in  a  constant 
electromagnetic  fieldo   VJe  find 

X  =  ut  , 

(3)      y  =  J  sin  (^yt  -  i  (  w  -  v')(  cos  ^t  -  1)      , 

z  =  ,^  (cos  6Jt  =  1)  +  i  (w  -  v'  )sin  Wt  +  v't  , 

where  tJ=  eB/mc 

The  time  of  flight  of  the  charged  particle  in  the  boundary  layer 
is  Immediately  found  to  be 

— -  arc  tan  ■,  —  , 
\U)\  V  =w 

where  the  i  sign  is  chosen  to  be  the  same  as  the  sign  of  the  charge  e 
on  the  particle  J  and  the  value  of  the  arc  tangent  is  chosen  between 
0  and  IT,   We  note  that  at  the  point  of  reentry  into  the  plasma,  the 
velocity  is  (u,-=Vpw)„   At  the  particle's  subsequent  emergence  into 
the  other  boundary  layer^  its  initial  velocity  (in  the  appropriate 
coordinate  system  there)  is  (u,Vj,-w)o   Consequently  the  time  of  flight 
there  is 

2       ^    +  V 

— •  arc  tan   ,  ■  ■  ■  » 

I  (x)  1         V  +  w 
Therefore  the  range  r  is  given  by 

(I4.)     r  =  -JL.  (arc  tan  -— -  +  arc  tan  -~-  )  . 
I  63  I         V  =w  V  +w 

After  this  result  is  Inserted  into  the  formula  (2)  for  the  loss 
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rate  Q,  it  is  just  a  matter  of  integration  to  determine  the  rate  at 
which  particles  are  lost  through  the  cusp„   In  general,  this  requires 
numerical  computation.   We  shall  be  able  to  evaluate  Q  explicitly  in 
a  fevj  special  cases. 

Of  course  the  answer  depends  on  the  value  of  E,  or,  equivalently, 
of  V' o   Ultimately  the  value  of  E  must  be  determined  by  an  analysis 
of  the  entire  space  charge  field  which  sxirroixnds  the  plasma  in  order 
to  equalize  the  rates  of  loss  of  ions  and  electronSo   In  the  remainder 
of  this  paper,  we  only  evaluate  the  loss  rate  Q,  for  certain  assigned 
values  of  Eo   We  note  that  the  loss  rates  for  electrons  and  ions  will 
turn  out  to  differ,  as  should  be  expected. 
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'-t-o   Loss  rs  zes   Q_j_ 

For  simplicity  we  now  take  a  distribution  of  particles  that  is 
raonoenergetic  as  well  as  isotroplco 

First  consider  the  case  E  =  Oo   The  range  r,  given  by  (Ll-), 
reduces  to  r  =  tcu/|6JI„   To  evaluate  the  loss  rate  Q,  which  is  given 
by  (2),  we  introduce  space  polar  coordinates; 

u  =  V  sin  0  cos  ^,      V  =  V  sin  9  sin  ^^      w  =  V  cos  9. 

The  angle  ■^    is  the  same  as  that  introduced  in  section  2j  see  Figiore  2b, 
The  speed  V  is  constants   The  distribution  function  f  becomes,  in 
these  variables, 

f(u,v,w)  du  dv  dw  =  Y—   sin  9  d9  d\l;, 

where  n  denotes  the  number  density^  Then  the  loss  rate  Q  is  found, 
by  a  simple  integration,  to  be 

(5)      Q   =   i  n  V  >.  ,  . 

where  \   =  V/jtOl  is  the  radius  of  gyration.   Electrons  and  ions  are 

2     2 

lost  at  the  same  rate  if  they  have  the  same  temperature,  n  V  =  m  V  . 

Although  the  electrons  move  faster,  the  size  of  the  "hole",  through 
which  they  leak  is  proportionately  smallerc  If  we  define  the  "hole 
size"  5  by  Q  =  n  u  5,  where  u  is  the  mean  speed 


u  =   I  u  f  du  dv  dw  /   I  f  du  dv  dw. 


u>  0  u  >0 

then  we  obtain  in  this  case 


5   =  ■»■  X 
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The  next  case  we  treat  is  more  important  in  that  it  provides  a 
partially  self-consistent  picture  of  the  boundary  layer.   We  adjust 
the  constant  value  of  the  electric  field  E  to  make  the  mean  penetra- 
tion of  ions  and  electrons  equal,   (The  Debye  length  is  assiomed  to 
be  small  compared  to  the  thickness  of  the  boundary  layer,) 

We  make  the  tentative  hypothesis  that  V   =  e/b  lies  midway 
between  V  and  V_,  i^e,.  that 

'/   «  v'  «  V   . 
+  — 

At  the  end  of  our  analysis,  we  shall  be  able  to  check  on  the  validity 
of  this  assumption.  For,  the  ions  the  range  r,  given  by  (Ij.),  in  this 
case  reduces  approximately  to 

(6j     r  =  -r--   arc  tan  —r   , 
^+         V 

Then,  using  the  trajectory  equations  (3)  it  is  easy  to  compute  d_^, 
the  mean  penetration  of  a  positive  ion  into  the  boundary  layer;  and 
we  obtain 

d   -1^ 

+   3  oJ  v' 
+ 

The  average  mean  penetration  d_j_s  averaged  over  the  positive  ion  dis- 
tribution, is  then  immediately  found  to  be 

I  -  -'" 

+  - 

On  the  other  hand,  for  electrons  we  ignore  V  =  E/b  in  comparison 
to  V  «  Then  we  find  the  mean  penetration  d_  of  an  electron  into  the 
boundary  layer  to  be  simply 
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/v  +  w 
d.  =  ,  ,  , — ■ 

The  average  mean  penetration  of  the  electrons  is  consequently 

d   -   -A. 

Now  we  set  d  =  d._  --  and  assuming  that  the  electron  and  ion 
temperatures  are  equal,  l„e<,   m_|_V ,  =  ^^^   '"   '"®  find 

V   =  ^  V   . 

For  deuteri-um,  we  approximately  have 

V;v';V  =i?r:l;7. 
+       -   10       '  ' 

confirming  our  hypothesis  about  the  relative  magnitudes  of  V^  V  , 
and  V_, 

This  analysis  results  in  a  boundary  layer  thickness  close  to  an 
electron  gyro-radius j,  and  the  value 

E  =  Ji.  V  B 
971  - 

for  the  electric  fieldo   We  are  now  able  to  apply  our  general  loss 
rate  formula  (2),   In  fact  since  E  hardly  alters  the  electron  orbitSj 
the  electron  loss  rate  is  given  as  before, 

(7)      Q„  =  3  n  V^  ?v_  , 

For  ions  we  use  for  the  range  r  ,  the  approximate  value 
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_  2u  V 
^+  v'  ' 

obtained  from  (6),   Inserting  this  expression  for  v_^   into  (2),  we 

find  for  the  ion  loss  rate 

^+  ~  32   V_  ''+  • 
The  "hole  size"  for  ions  is  then 

whereas  for  electrons  it  is 

Therefore  we  have 

5^  I  5_   =  2.7  . 

For  deuterium,  we  obtain 

Q^s  Q_  =   1  :  35  . 

The  difference  in  loss  rates  must  be  equalized  by  external  space 
charge  effects.   For  a  qualitative  resolution  of  this  difficulty, 
see  [I]-]  .   The  equalized  loss  rate  is  intimately  connected  with  the 
thickness  of  the  boundary  layer,  and  there  is  a  factor  of  at  least 
35  at  stake.   In  the  picture  presented  here,  the  electron  orbits  are 
hardly  affected  by  the  electric  field  in  the  boxondary  layer,  and  the 
ion  orbits  are  pulled  in  almost  to  the  electron  orbits. 
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